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Abstract 



A realistic microscopic calculation of the unpolarized quark Generalized 
Parton Distribution (GPD) of the "^He nucleus is presented. In Impulse 
Approximation, is obtained as a convolution between the GPD of the 
internal nucleon and the non-diagonal spectral function, describing properly 
Fermi motion and binding effects. The proposed scheme is valid at low values 
of A^, the momentum transfer to the target, the most relevant kinematical 
region for the coherent channel of hard exclusive processes. The obtained 
formula has the correct forward limit, corresponding to the standard deep 
inelastic nuclear parton distributions, and first moment, giving the charge 
form factor of ^He. Nuclear effects, evaluated by a modern realistic potential, 
are found to be larger than in the forward case. In particular, they increase 
with increasing the momentum transfer when the asymmetry of the process is 
kept fixed, and they increase with the asymmetry at fixed momentum transfer. 
Another relevant feature of the obtained results is that the nuclear GPD 
cannot be factorized into a A^-dependent and a A-^-independent term, as 
suggested in prescriptions proposed for finite nuclei. The size of nuclear effects 
reaches 8 % even in the most important part of the kinematical range under 
scrutiny. The relevance of the obtained results to study the feasibility of 
experiments is addressed. 
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I. INTRODUCTION 



Generalized Parton Distributions (GPDs) [1] parametrize the non-perturbative hadron 
structure in hard exclusive processes (for comprehensive reviews, see, e.g., [2-5]). The 
measurement of GPDs would provide information which is usually encoded in both the elastic 
form factors and the usual Parton Distribution Functions (PDFs) and, at the same time, it 
would represent a unique way to access several crucial features of the miclcon [6,7], such as 
its angular momentum content [7] and its structure in the transverse plane [8]. According 
to a factorization theorem derived in QCD [9], GPDs enter the long-distance dominated 
part of exclusive lepton Deep Inelastic Scattering (DIS) off hadrons. In particular. Deeply 
Virtual Compton Scattering (DVCS), i.e. the process eH — > e'H'^ when ^ rrijj, is one 
of the the most promising to access GPDs (here and in the following, is the momentum 
transfer between the leptons e and e', and the one between the hadrons H and H') 
[6,7,10]. Therefore, relevant experimental efforts to measure GPDs by means of DVCS off 
hadrons are hkely to take place in the next few years. As a matter of fact, a few DVCS data 
have been already published [11,12]. 

Recently, the issue of measuring GPDs for nuclei has been addressed. The first paper on 
this subject [13], concerning the deuteron, contained already the crucial observation that the 
knowledge of GPDs would permit the investigation of the short light-hke distance structure 
of nuclei, and thus the interplay of nucleon and parton degrees of freedom in the nuclear 
wave function. In standard DIS off a nucleus with four-momentum Pa and A nucleons of 
mass M, this information can be accessed in the region where Axbj — 2Mj^ > being 
XBj = Q'^H'^Pa ■ q) and the energy transfer in the laboratory system. In this region 
measurements are very difficult, because of vanishing cross-sections. As explained in [13], 
the same physics can be accessed in DVCS at much lower values of Xsj- The usefulness of 
nuclear GPDs has been stressed also for finite nuclei in Ref. [14], where it has been shown 
that they could provide us with peculiar information about the spatial distribution of energy, 
momentum and forces experienced by quarks and gluons inside hadrons. Since then, DVCS 
has been extensively discussed for nuclear targets. Impulse Approximation (lA) calculations, 
supposed to give the bulk of nuclear effects at 0.05 < Axbj < 0.7, have been performed for 
the deuteron [15] and for spinlcss nuclei [16]. For nuclei of any spin, estimates of GPDs have 
been provided and prescriptions for nuclear effects have been proposed in [17]. Recently, an 
analysis of nuclear DVCS has been performed beyond lA, with estimates of shadowing effects 
and involving therefore large light-like distances and correlations in nuclei [18]. It was found 
that these effects are sizable up to Axbj — 0.1, i.e. up to larger values of Axbj with respect 
to normal DIS. Besides, the possibility of measuring DVCS at an electron-ion-collider has 
also been established in [18]. 

The study of GPDs for ^He is interesting for many aspects. In fact, ^He is a well known 
nucleus, for which realistic studies are possible, so that conventional nuclear effects can be 
safely calculated. Strong deviations from the predicted behaviour could be therefore as- 
cribed to exotic effects, such as the ones of non-nucleonic degrees of freedom, not included 
in a realistic wave function. Besides, ^He is extensively used as an effective neutron target. 
In fact, the properties of the free neutron are being investigated through experiments with 
nuclei, whose data are analyzed taking nuclear effects properly into account. Recently, it has 
been shown that unpolarized DIS off three body systems can provide relevant information on 
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PDFs at large xsj [19-22], while it is known since a long time that its particular spin struc- 
ture suggests the use of "^/Je as an effective polarized neutron target [23-26]. Polarized ^He 
will be therefore the first candidate for experiments aimed at the study of spin- dependent 
GPDs in the free neutron, to unveil details of its angular momentum content. 

In this paper, an Impulse Approximation (lA) calculation of the quark unpolarized GPD 
of "^He is presented. A convolution formula is derived and afterwards numerically eval- 
uated using a realistic non-diagonal spectral function, so that Fermi motion and binding 
effects are rigorously estimated. The proposed scheme is valid for -C Q^, and despite 
of this it permits to calculate GPDs in the kinematical range relevant to the coherent, no 
break-up channel of deep exclusive processes off ^He. In fact, the latter channel is the most 
interesting one for its theoretical implications, but it can be hardly observed at large A^, 
due to the vanishing cross section [18]. The nuclear GPDs obtained here are therefore a 
prerequisite for any calculation of observables in coherent DVCS off ^He, although they can 
be compared neither with existing data, nor with forth-coming ones. A detailed analysis of 
DVCS off ^He, with estimates of observables, such as cross-sections or spin asymmetries, is 
in progress and will follow in a separate publication. Thus, the main result of this investiga- 
tion is not the size and shape of the obtained for ^He, but the size and nature of nuclear 
effects on it. This will permit to test directly, for the ^He target at least, the accuracy of 
prescriptions which have been proposed to estimate nuclear GPDs [17], providing a useful 
tool for the planning of future experiments and for their correct interpretation. 

The paper is organized as follows: The notation used is introduced in section 2, together 
with the derivation of the main result, a convolution formula obtained from an lA analysis. 
Such a formula is used to evaluate the numerical results shown in section 3 and thoroughly 
discussed in section 4, where the most interesting outcome of the present study is to be 
found. Eventually, conclusions are drawn in the fifth section. 



The formalism introduced in Ref. [2] is adopted. One has to think to a spin 1/2 hadron 
target, with initial (final) momentum and helicity P{P') and s(s'), respectively. The GPDs 
ifq(a;, ^, A^) and Eq[x,^, A^) are defined through the expression 



where A = P' — P is the 4-momentum transfer to the hadron, is the quark field and M 
is the hadron mass. It is convenient to work in a system of coordinates where the photon 
4-momentum, = (go, q), and P = {P + P')/2 are coUinear along z. The ^ variable in the 
arguments of the GPDs is the so called "skewedness" , parametrizing the asymmetry of the 
process. It is defined by the relation 



II. FORMALISM 




U{P,s) , (1) 
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where n is a light-like 4- vector satisfying the condition n-P = 1. One should notice that the 
variable ^ is completely fixed by the external lepton kinematics. Here and in the following, 
use is made of the definition = (a° ± a^)/\/2. As explained in [7], GPDs describe the 
amplitude for finding a quark with momentum fraction x + ^ (in the Infinite Momentum 
Frame) in a hadron with momentum (1 + and replacing it back into the nucleon with a 
momentum transfer A. Besides, when the quark longitudinal momentum fraction x of the 
average nucleon momentum P is less than —C,, GPDs describe antiquarks; when it is larger 
than ^, they describe quarks; when it is between and ^, they describe qq pairs. One 
should keep in mind that, in addition to the variables x, ^ and exphcitly shown, GPDs 
depend, as the standard PDFs, on the momentum scale at which they are measured or 
calculated. Such a dependence will not be discussed in the present paper and, for an easy 
presentation, it will be omitted. The values of ^ which are possible for a given value of A^ 
are: 

< ^ < v^=A2/V4M2 - A2 . (3) 

The well known natural constraints of Hq{x,^, A'^) are: 

i) the so called "forward" or "diagonal" hmit, P' — P, i.e., A^ = ^ = 0, where one 
recovers the usual PDFs 

Hg{x,0,0)=q{x) ; (4) 

ii) the integration over x, yielding the contribution of the quark of fiavour q to the Dirac 
form factor (f.f.) of the target: 

JdxH,ix,i,A')=F!{A') ; (5) 

iii) the polynomiality property [2], involving higher moments of GPDs, according to 
which the integrals of x'^W and of x'^E^ are polynomials in ^ of order n + 1. 

In [27,28], an expression for Hq{x, ^, A^) of a given hadron target, for small values of 
has been obtained from the definition Eq. (1). 
It reads: 

i/,(x,e,A2)= Jdksl^'^-ix + O^ X 

(6) 



X 



T.iP'lbUk'- + A\ + A^) k^)\P) 



{2nfVko ^ 

where states and creation and annihilation operators are normalized as follows 

(P'|P) = (27r)^5(P'+ - P+)(5(P1 - Pa) (7) 

and 

A^x), h\k'+, fcl)} = {2nf6{k'+ - k+)6{k'^ - k^) , (8) 

respectively. 
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All the steps leading from Eq. (1) to Eq. (6) are fully described in Ref. [28]. They 
include initially the use of light-cone spinors and fields, neglecting terms of order in 
the right hand side of Eq. (1), and later a transition to the expression (6), where the same 
quantities are normalized according to Eqs. (7) and (8). 

The case of a spin 1/2 hadron with three composite constituents has been discussed in 
[28]. In that paper, the GPDs of the proton have been studied, assuming that it is made of 
complex constituent quarks. In here, the approach will be extended to ^He. The GPD 
of ^He will be obtained in lA as a convolution between the non-diagonal spectral function 
of the internal nucleons, and the GPD of the nucleons themselves. 

The scenario is depicted in Fig. 1 for the special case of coherent DVCS, in the handbag 
approximation. One parton with momentum A;, belonging to a given nucleon of momentum 
p in the nucleus, interacts with the probe and it is afterwards reabsorbed, with momentum 
A; -|- A, by the same nucleon, without further re-scattering with the recoiling system of 
momentum Pr. Finally, the interacting nucleon with momentum p -|- A is reabsorbed back 
into the nucleus. The analysis suggested here is quite similar to the usual I A approach to 
DIS off nuclei [30,31]. 

In the class of frames discussed above, and in addition to the kinematical variables, x 
and ^, already defined, one needs a few more to describe the process. In particular, x' and 

for the "internal" target, i.e., the nucleon, have to be introduced. The latter quantities 
can be obtained defining the components of the momentum k and /c -|- A of the struck 
parton before and after the interaction, with respect to P+ and p+ = |(p + 

k+^{x + i)P+^{x' + OP^ , (9) 
{k + A)+ = (x - = {x' - . (10) 

From the above expressions, ^' and x' are immediately obtained as 

x' = ^x (12) 
and, since ^ = — A+/(2P+), if z =p+/P+, one also has 

These expressions have been already found and used in [13,16,17,28]. 

In order to derive a convolution formula in lA for H^, the procedure used for standard 
DIS will be adopted [30,31]. First, in Eq. (6), two complete sets of states, corresponding to 
the interacting nucleon and to the recoiling two-body system, are properly inserted to the 
left and right-hand sides of the quark operator: 

Hl{x,^,A') = {P'S\ ^ {\PkS'n)\p's')}{{PkS'R\¥s'\} 

^^^^^^ ' ^ E^Ia(^^ + a^^± + a^)M^^^±J 



I 



(27r)3 k 




E {\PrSr)\ps)}{{PrSr\{ps\} \PS); (14) 

Pr,Sr,p,s 
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second, since eventually use has to be made of NR nuclear wave functions, a state |p) has 
to be normalized in a NR manner: 

im-c^^w-p}, (15) 

so that in Eq. (14) one has to perform the substitution: 

\P>^\l^\P> ■ (16) 
Since, using lA in the intrinsic frame of ^He one has, for the NR states: 

{{PrSr\{ps\}\PS) = {PRSR,ps\PS){2nfd\P - Pr- p)ds,sus , 
a convolution formula is readily obtained from Eq. (14): 



Hlix, A')c^j:[dE[ dp Pl^^±^Pl,(p,p + A, E)^-H^{x', C', A^) + O (e) = (17) 

N \ Po[P + ^jo ? 

= T.JdEl dp[PUp,p + A, E) + 0{f/M\ AVM^)] 



N 



X ^H^{x',^',A') + 0(e) ■ (18) 



In the above equations, the kinetic energies of the residual nuclear system and of the recoiling 
nucleus have been neglected, and -P^(p,p + A, £■) is the one-body off-diagonal spectral 
function for the nucleon N in ^He: 

PUp,P + a, E) = ^Wj2{P'M\{P - p)Sr, {p + A)s){{P - p)Sr,ps\PM) x 

X 5{E - E^i^ - E*j,) . (19) 

Besides, the quantity 



//f(x',r,A^) = ldks(^^-{x'+e)j 



{27rrVko , 

ydks(^^-{x+o) 



(2^|^E(p'l^lA(^+ + A+,^x + A^)^,A(A;+,Mb) (20) 

is, according to Eq. (6), the GPD of the bound nucleon N up to terms of order 0{^^), and 
in the above equation use has been made of Eqs. (11) and (12). 

The delta function in Eq (18) defines E, the so called removal energy, in terms of E^nin — 
\E3Hf,\ — \E2h\ = 5.5 MeV and E^, the excitation energy of the two-body recoiling system. 
The main quantity appearing in the definition Eq. (19) is the overlap integral 



6 



{PM\PrSrJs) = j dye'^^-y{r,^'jf{x)\^^{x,^) , (21) 

between the eigenfunction of the ground state of ^He, with eigenvalue Esh^ and third 
component of the total angular momentum M, and the eigenfunction "^r^, with eigenvalue 
Er — E2 + E*j^ of the state R of the intrinsic Hamiltonian pertaining to the system of two 
interacting nucleons [32] . Since the set of the states R also includes continuum states of the 
recoiling system, the summation over it! involves the deuteron channel and the integral over 
the continuum states. 

Concerning Eqs (18-21), two comments are in order. 

The first concerns the accuracy of the actual calculations which will be presented. In the 
following, a NR spectral function will be used to evaluate Eq. (18), so that the accuracy of 
the calculation is of order O In fact, if use is made of NR wave functions 

to calculate the non-diagonal spectral function, the result holds for p^, [p + A)^ << M^. 
The same constraint can be written p^, A^ << M^. While the first of these conditions is the 
usual one for the NR treatment of nuclei, the second forces one to use Eq. (18) only at low 
values of A^, for which the accuracy is good enough. As it will be explained, the interest of 
the present calculation is indeed to investigate nuclear effects at low values of A^. 

Second, from the study of forward DIS off nuclei, it is known that, to properly estimate 
nuclear effects, in going from a covariant formalism to one where use can be made of the usual 
nuclear wave functions, one has to keep the correct normalization [30,33]. This procedure 
leads to the appearance of the "fiux factor", represented in Eq. (17) by the expression 

\lp^{p + A)+/po(p + A)o (which gives, in forward DIS, the usual p'^/po term). This factor 

gives one at order O (jP/M'^, A^/M^^, which is the accuracy of the present analysis, and it 
will not be included in the actual calculation. In other words, if one takes into account that 
the intrinsic light-front frame and the usual rest frame of ^He, where the wave functions 
are evaluated, are not the same frame, the flux factor should be introduced [30,33]; anyway, 
in the NR approach used here, the two frames do not differ. This is an effect of the used 
approximations, which have to be relaxed if one wants to be predictive in more general 
processes at higher momentum transfer. In fact, in that case the accuracy of Eq. (18) is not 
good enough any more and the calculation, in order to be consistent, has to be performed 
taking into account relativistic corrections. Work is being done presently in this direction, 
and the discussion of this important point will be included in a following paper, describing 
a calculation performed in a light-front framework. One final remark about the fiux-factor: 
to neglect the corresponding quantity for the quark inside the internal nucleon, i.e. the term 
k'^/ko in the GPD of the bound nucleon, Eq. (20), would be a rough approximation, with 
respect to the nuclear case, because the motion of quarks in the nucleon is relativistic, even 
at the constituent level. In any case in the present calculation, as described in the following 
section, the GPD for the bound quark will not be evaluated, while an available model for 
it will be used. The main emphasis of the present approach, as already said, is not on the 
absolute values of the results, but in the nuclear effects, which can be estimated by taking 
any reasonable form for the internal GPD. 
Eq. (18) can be written in the form 

H'^{x,^,A')=J2 [ dE [ dpPUp,p + A)^H^{x',e,^') = 
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-^CiJ<iElm,P + A)6[.-i)H^[^.A'^) . (22) 

Taking into account that 

^- I = [5(1 + = ^ + 1^(1 + = ^ + (23) 
Eq. (22) can also be written in the form: 

H'^ix, A^) = E /" ^hUz, f , A^) , (24) 



where 



(25) 



One should notice that Eqs. (24) and (25) or, which is the same, Eq. (18), fulfill the 
constraint i) — in) previously listed. 

The constraint i), i.e. the forward limit of GPDs, is certainly verified. In fact, by 
taking the forward limit (A^ ^ 0,.^ — 0) of Eq. (24), one gets the expression which is 
usually found, for the parton distribution 53(2;), in the lA analysis of unpolarized DIS off 
^He [20,30,31]: 

qs{x) = H'^ix, 0, 0) = E /' -fUz) qN (-) ■ (26) 

Jx Z \Z / 

In the latter equation, 

fUz) = h%{z, 0, 0) ^JdEj dpPUp, E)6 (^z - ^ j (27) 

is the light-cone momentum distribution of the nucleon N in the nucleus, qN{x) = H^[x, 0, 0) 
is the distribution of the quark of fiavour q in the nucleon and P^{p, E), the A^ — > 
hmit of Eq. (22), is the one body spectral function. 

The constraint ii), i.e. the x— integral of the GPD Hq, is also naturally fulfilled. In fact, 
by x— integrating Eq. (24), one easily obtains: 

= E / dx'H^ix',e,A') [ dzh%{z,^,A') = 
^Y:F^{A')F%{A')^FI{A'). (28) 

N 

In the equation above, (A^) is the contribution, of the quark of fiavour g, to the nuclear 
f.f.; Fj^(A^) is the contribution, of the quark of fiavour g, to the nucleon N f.f.; F^(A^) is 



the so-called ^He "pointlike f.f.", which would represent the contribution of the nucleon N 
to the f.f. of ^He if were point- like. F^(A^) is given, in the present approximation, by 

F3(A2) = jdEj dpPUp,P + A,E) = I dzhUz,^,A') . (29) 

Eventually the polynomiality, condition in), is formally fulfilled by Eq. (18), although 
one should always remember that it is a result of order C(^^), so that high moments cannot 
be really checked. 

Summarizing this section, we have derived in lA a convolution formula for the GPD 
of ^He, Eq. (18) (or, which is equivalent, Eq. (24)), at order O (^-^, in terms of 

a non-diagonal nuclear one-body spectral function, Eq. (19), and of the GPD of the 
internal nucleon. 



III. NUMERICAL RESULTS 

H^{x,^j A^), Eq. (18), has been evaluated in the nuclear Brcit Frame. 

The non-diagonal spectral function Eq. (19), appearing in Eq. (18), has been calculated 
along the hues of Ref. [34], by means of the overlap Eq. (21), which exactly includes the 
final state interactions in the two nucleon recoiling system, the only plane wave being that 
describing the relative motion between the knockcd-out nucleon and the two-body system 
[32]. The realistic wave functions and "^'^ in Eq. (21) have been evaluated using the 
AV18 interaction [35] and taking into account the Coulomb repulsion of protons in ^He. In 
particular has been developed along the hues of Ref. [36]. The same overlaps have been 
already used in Ref. [20]. 

The other ingredient in Eq. (18), i.e. the nucleon GPD , has been modelled in 
agreement with the Double Distribution representation [29], as described in [37]. For the 
reader convenience, the explicit form of is listed below [3,29]: 

H^{x, t A^) = dx S{x + ^a- x)%{x, a, A'^)da . (30) 

J-l J-l+\x\ 

With some care, the expression above can be integrated over x and the result is explicitly 
given in [3]. 

In [29], a factorized ansatz is suggested for the DD's: 

^,{x, a, A^) = h,{x, a, A^)%{x)Fg{A^) , (31) 

with the a dependent term, hg{x,a, A^), which has the character of a mesonic amplitude, 
fulfilling the relation: 

fl-\x\ 

/ hg{x,a,A'^)da = l . (32) 

J-l+\x\ 

Besides, in Eq. (31) ^q{x) represents the forward density and, eventually, Fg(A^) the con- 
tribution of the quark of flavour q to the nucleon form factor. Eq. (30) fulfills the crucial 
constraints of GPDs, i.e., the forward limit, the first-moment and the polynomiality condi- 
tion. One needs now to model the three functions appearing in Eq. (31). For the amplitude 
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hq, use will be made of one of the simple normalized forms suggested in [29], on the bases 
of the symmetry properties of DD's: 

For the forward distribution ^g{x), the general simple form 

= (34) 

with a = 0.5, has been taken. 

Concerning the model used for the internal nucleon , two caveats are in order. 

First of all, Eq. (30), with the choices Eqs. (33) and (34), corresponds to the NS, valence 
quark contribution proposed in [37], which is symmetric in the variable x. 1 stress again 
that the main point of the present study is not to produce realistic estimates for observables, 
but to investigate and discuss nuclear effects, which do not depend on the form of any well- 
behaved internal GPD, whose general structure is safely simulated by Eqs. (30) - (34). 
For the moment being, only the NS part of is therefore calculated; if one wanted to 
estimate the DVCS cross-sections, also the quark singlet and the gluon contributions would 
be certainly necessary. As already said, a detailed study of DVCS on ^He is in progress and 
will be shown elsewhere. 

Secondly, to take advantage of the used parametrization, Eqs. (30) - (34), of the nucle- 
onic GPD, one should work in a symmetric frame for the nucleon, such as the Breit Frame. 
Actually, the present calculation is performed in the nuclear Breit Frame, which does not 
coincide with the Breit Frame of the internal nucleon. For this reason, one should perform a 
boost to properly work in the nuclear Breit Frame. Anyway, as it is pointed out in Ref. [16], 
any NR evaluation of nuclear GPDs, as a result of the NR reduction, is intrinsically frame 
dependent and it is valid only with accuracy 0{p'^/M^), like Eq. (18). In this calculation 
therefore, the effect of using a different reference frame for the nucleon is not taken into 
account, since in any case it would be an effect of order 0(p^/M^), which is the overall 
accuracy of the present work [cf. Eq. (18)]. In the above procedure, therefore, no further 
approximations are introduced in addition to the ones leading to Eq. (18). To overcome 
such a problem, one could use a light-front approach [38], as it is done for the GPDs of the 
deuteron in [13,15]. 

Eventually, the -^^(A^) term in Eq. (31), i.e. the contribution of the quark of flavour q 
to the nucleon form factor, has been obtained from the experimental values of the proton, 
Fi, and of the neutron, F^, Dirac form factors. For the u and d flavours, neglecting the 
effect of the strange quarks, one has 

F„(A2) = i(2Ff(A^) + Ff(A^)), 

Fd{A^) = 2F{\/\^) + Ff (A^) . (35) 
The contributions of the flavours u and d to the proton and neutron f.f. are therefore 

F„^(A^) = ^F„(A^) , 

F| = -^F,(A2) (36) 
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and 

F2{A') = —F^iA') , (37) 

respectively. 

For the numerical calculations, use has been made of the parametrization of the nucleon 
Dirac f .f . given in Ref . [39] . 

Now the ingredients of the calculation have been completely described, so that numerical 
results can be presented. First of all, the forward limit of H^, Eq. (26), will be discussed, 
together with its a;— integral, Eq. (28). 

In Fig. 2, it is shown the forward limit of the ratio 

H^(x,^,A'^) 

^' ^ 2i/f(a;,e,A2)Vi/«(a;,e,A2) ' ^^^^ 

i.e. the quantity 

R(.OO) ""^'^"'"'"^ ^'("^ (39) 

' ' ^ 2H^(x,0,0) + H^(x,0,0) 2qP(x)+q-(x) ^ ^ 

for the flavour q — u,d, as a function of — 3x. In the above equation, the numerator is 
given by Eq. (26), while the denominator clearly represents the distribution of the quarks 
of flavour q in ^He if nuclear effects are completely disregarded, i.e., the interacting quarks 
are assumed to belong to free nucleons at rest. The behaviour which is found is typically 
i^MC— like, as it is usually obtained in lA studies of DIS on nuclei [20,30,31], so that, in the 
forward limit, well-known results are recovered. One should notice that, had the structure 
functions ratio, instead of the PDFs one, been shown, the ratio would be 1 for x = 0, due to 
the normalization of the spectral function in Eq. (27). It is also useful, for later convenience, 
to realize that nuclear effects for the d flavour are a bit larger than those which are found 
for the u flavour. This is due to the fact that the forward d distribution is more sensitive 
than the u one to the neutron light-cone distribution, Eq. (27), which is different from that 
of the proton. In fact, the average momentum of the neutron in ^He is larger than the one 
of the proton [32,31]. 
In Fig. 3, the quantity 

J dxHl^x^i^') = \'£'£jdxj ^h%{z,C,A^)H^ (l'!'^') ' 

is shown. According to Eq. (28), in the present approximation, this should give, up to con- 
tributions of heavier quarks, the charge form factor Fp^^( A^) of ^He (the usual normalization 
F4(0) = 1 is chosen): 

F!,{A') = 1[F^{A') + F!{A')]. (41) 

In the figure, the calculated integral is compared with experimental data of |Fp^^(A^)| [40]. 
It is found that the present approach reproduces well the data up to a momentum transfer 
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— = 0.25 GcV^, which is enough for the aim of this calculation. In fact, the region of 
higher momentum transfer is not considered here, being phenomenologically not relevant for 
the calculation of GPDs entering coherent DVCS. The full curve in the left panel of Fig. 3 
is very similar to the one-body contribution to the ^He charge f.f. shown in [41], as it must 
be, due to the lA used here. The agreement with data in the relevant kinematical region, 
— < 0.25 GeV^, confirms that the inclusion of two-body currents is not required in the 
present calculation. 

As an illustration, the result of the evaluation of H^{x, ^, A^) and H^{x, ^, A^) by means 
of Eq. (18) is shown in Figs. 4 and 5, for A^ = -0.15 GeV^ and A^ = -0.25 GeV^ 
respectively, as a function of Xs and ^3 = 3^. The GPDs are shown for the ^3 range allowed 
by Eq. (3) and in the ^3 > region. The X3 < one, being symmetric to the latter, is not 
interesting and not shown. 

The quality and size of the nuclear effects are discussed in the next section. 



IV. DISCUSSION OF NUCLEAR EFFECTS 

The full result for the GPD H^, Eq. (18), shown in Figs. 4 and 5, will be now compared 
with a prescription based on the assumptions that nuclear effects are completely neglected 
and the global A^ dependence can be described by the f.f. of ^He: 

H}^'\x,^,A') = 2H}^{x,^,A') + H}-{x,^,A') , (42) 

where the quantity 

H}^(x, A^) = H^(x, OF^i^') (43) 

represents the flavor q effective GPD of the bound nucleon N = n,p in ^He. Its x and ^ 
dependences, given by the function [x,^), is the same of the GPD of the free nucleon N 
(represented in this calculation by Eq. (30)), while its A^ dependence is governed by the 
contribution of the quark of flavor q to the ^He f.f., F^{A?'). 

The effect of Fermi motion and binding can be shown through the ratio 

i.e. the ratio of the full result, Eq. (18), to the approximation Eq. (42). The latter is 
evaluated by means of the nucleon GPDs used as input in the calculation, and taking 

^.^(A^) = f ^i(A^) , (45) 



F|(A^) = -\fI{A^) , (46) 

where Fjj(A^) is the f.f. which is calculated within the present approach, by means of Eq. 
(28). The coefficients 10/3 and —4/3 are simply chosen assuming that the contribution of 
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the valence quarks of a given flavour to the f.f. of ^He is proportional to their charge. One 
should remember that the normalization of the f.f. has been chosen in Eq. (41). 

The choice of calculating the ratio Eq. (44) to show nuclear efi'ects is a very natural 
one. As a matter of fact, the forward limit of the ratio Eq. (44) is the same of the ratio 
Eq. (38), yielding the EMC-hke ratio for the parton distribution q and, if ^He were made 
of free nuclcon at rest, the ratio Eq. (44) would be one. This latter fact can be immediately 
realized by observing that the prescription Eq. (42) is exactly obtained by placing z = 1, 

1. e. no Fermi motion effects and no convolution, into Eq. (18). In the present situation 
the ratio Eq. (38) cannot be used to show nuclear effects, since it does not contain the 
dependence coming from the nuclear structure, i.e. the one provided by the non-diagonal 
spectral function in the present approach [cf. Eq. (18)]. 

One should note that the prescription suggested in Ref. [17] for finite nuclei in the valence 
quark sector, basically assuming that the nucleus is a system of almost free nucleons with 
approximately the same momenta, has the same dependence of the prescription Eq. (42). 

Results are presented in Fig. 6 and 7, where the ratio Eq. (44) is shown for A^ = —0.15 
GeV^ and A^ = —0.25 GeV^, respectively, as a function of X3, for three different values of 
^3, for the fiavours u and d. 

Some general trends of the results are apparent: 

i) nuclear effects, for < 0.7, are as large as 15 % at most. Larger effects at higher 
values are actually related with the vanishing denominator in Eq. (44). 

ii) Fermi motion and binding have their main effect for x^ < 0.3, at variance with what 
happens in the forward limit (cf. Fig. 2). 

iii) at fixed A^, nuclear effects increase with increasing for X3 < 0.3. 

iv) at fixed ^, nuclear effects increase with increasing A^, for 2:3 < 0.3. 

v) nuclear effects for the d flavour are larger than for the u flavour. 

The behaviour described above can be explained as follows. As already said in section 

2, in lA and in the forward limit, at 0:3 = one basically recovers the spectral function 
normalization and no nuclear effects, so that the ratio Eq. (38) slightly differs from one. 
This is not true of course in the present case, due to nuclear effects hidden not only in 
the x' dependence, but also in the ^' one, according to its deflnition, Eq. (13). Moreover, 
even if X3 = ^3 = 0, in the present situation the ratio Eq. (44) does not give the spectral 
function normalization as in the forward case, because of the A^ dependence. One source 
of such dependence is that, in the approximation Eq. (42), it is assumed, through Eqs. 
(45) and (46) that the quarks u and d, belonging to the protons or to the neutron in ^ife, 
contribute to the charge f.f. in the same way, being the contribution proportional to their 
charge only. Actually, the effect of Fermi motion and binding is stronger for the quarks 
belonging to the neutron, having the latter a larger average momentum with respect to the 
proton [32,31]. This can be seen noticing that the pointlike f.f., Eq. (29), for the proton, 
shows a stronger A^-dependence, with respect to the neutron one, the difference being 17 
% (23 %) at A^ = -0.15 GeV^ (A^ = -0.25 GeV^). The prescription given by Eqs. (45) 
and (46) could be correct only if the pointlike f.f. had a similar A^ dependence. Besides, 
nuclear effects studied by means of the ratio Eq. (44) at fixed x and depend on A^, 
showing clearly that such a dependence cannot be factorized, i.e. the nuclear GPD cannot 
be written as the product of a A^ dependent and a A^ independent term, confirming what 
has been found for the deuteron case in Ref. [15]. One should notice that, if factorization 
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were valid, Figs. 6 and 7 would be equal. This fact clearly indicates that a model based on 
the assumption of factorization, such as the one of Ref. [17], is not motivated and cannot be 
used to parametrize nuclear GPDs for estimates of DVCS cross sections and asymmetries 
for light nuclei. 

The fact that nuclear effects are larger for the d distribution is also easily explained 
in terms of the different contribution of the spectral functions for the protons and the 
neutron, the latter being more important for the GPDs of the d rather than for the ones 
of the u flavour. This has been already shown in the forward case in Fig. 2, and shortly 
commented in the previous section. When the A^-dependence is studied, the effect found in 
the forward case increases, being impossible, as explained above, to really obtain the u and 
d contributions to the ^He f.f. by means of a simple charge rescaling, as it is done in Eqs. 
(45) and (46). In particular, the approximation Eq. (46) for the d flavour is worse than the 
approximation Eq. (45) for the u flavour, the disagreement increasing with increasing A^. 
One finds that F^(A^) given by Eq. (45) differs from the one calculated through Eq. (28) 
by 4 % (6 %) at A^ = -0.15 GeV^ (-0.25 GeV^), while for the d flavour the disagreement 
reaches 9 % (13 %) for the same value of the momentum transfer. This flavour dependence 
is to be expected for any target with isospin different from zero. 

As already said, the calculation of DVCS observables is beyond the scope of the present 
paper and will be presented elsewhere. Anyway, a first rough estimate of nuclear effects on 
DVCS observables can be already sketched. In fact it is known that the point x — ^ gives 
the bulk of the contribution to hard exclusive processes, since at leading order in QCD the 
amplitude for DVCS and for meson electroproduction just involve GPDs at this point, which 
is therefore also the easiest region to access experimentally. In order to figure out how Fermi 
motion and binding affect the "slice" H^{^, ^, A^), the ratio (44) for x = ^, is shown in Fig. 
8, as a function of ^3 = ^3, for the flavor u. It can be seen that even in this crucial region 
nuclear effects arc systematically underestimated by the approximation Eq. (42), up to a 
maximum of 8% for the flavor u. 

The issue of applying the obtained GPDs to calculate DVCS off ^'He, to estimate cross- 
sections and to establish the feasibility of experiments, as it has been done already for the 
deuteron target in Ref. [15], is in progress and will be presented elsewhere. In particular, 
the contributions of the incoherent break-up channel, and of the shadowing effects, beyond 
lA, at ^ ~ XBj/'i < 0.05, have to be evaluated. Besides, the study of polarized GPDs will 
be very interesting, due to the peculiar spin structure of ^He and its implications for the 
study of the angular momentum of the free neutron. 

V. CONCLUSIONS 

In this paper, a realistic microscopic calculation of the unpolarized quark GPD 
ifg(a;, ^,A^) for ^He is presented. In an Impulse Approximation framework, a convolu- 
tion formula has been obtained where Fermi motion and binding effects are properly taken 
into account through an off-diagonal spectral function. The range of validity of the scheme 
is |A^| < 0.25 GeV^, the most relevant for the coherent channel of hard exclusive processes. 
The forward hmit and the x-integral of the obtained H^{x, ^, A^) are formally and numeri- 
cally in agreement with the theoretical and experimental knowledge of the ^He nucleus. 
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Nuclear effects are found to be larger than in the forward case and to increase with 
at fixed ^, and with ^ at fixed A^. In particular the latter A^ dependence does not simply 
factorize, in agreement with previous findings for the deuteron target and at variance with 
prescriptions proposed for finite nuclei. Moreover, nuclear effects are found to be larger for 
the d flavour than for the u one, being the target non isoscalar. In particular, it has been 
shown that nuclear effects are as large as 8 %, even at the crucial point x = ^. 

The shown results represent a prerequisite for any calculation of DVCS observables. 
The evaluation of these quantities, to establish the feasibility of experiments, is in progress 
and will be presented in a later paper, together with the contributions of the incoherent, 
break-up channel, and of shadowing effects, beyond lA, at small values of ^ ~ a^sj/S. The 
present analysis is also a first step towards a complete study of the GPDs for ^He, which 
could provide, in the polarized case, with an important tool to unveil details of the angular 
momentum structure of the free neutron. 
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FIGURE CAPTIONS 



Fig. 1: The handbag contribution to the coherent DVCS process off ^He, in the present 
approach. 

Fig. 2: The ratio Eq. (39) as a function of xs, for the flavour u (left panel) and d (right 
panel) . 

Fig. 3: The charge f.f. of ^He^ calculated by means of Eq (41) (full line), is compared with 
a parametrization of data [40] (dashed line), in the region of low-A^, relevant to the present 
study. 

Fig. 4: In the left panel, for the ,^3 values which are allowed at = —0.15 GeV^ according 
to Eq. (3), Hl{x^,i^,/^^), evaluated using Eq. (24), is shown for 0.05 < 0:3 < 0.8. The 
symmetric part at ^3 < is not presented. In the right panel, the same is shown, for the 
flavour d. 

Fig. 5: The same as in Fig. 4, at = -0.25 GeV^. 

Fig. 6: In the left panel, the ratio Eq. (44) is shown, for the u flavour and A^ = —0.15 

GeV^, as a function of x^. The full line has been calculated for ,^3 = 0, the dashed line 
for ,^3 = 0.1 and the long-dashed one for ,^3 = 0.2. The symmetric part at 0:3 < is not 
presented. In the right panel, the same is shown, for the flavour d 

Fig. 7: The same as in Fig. 6, at A^ = -0.25 GeV^. 

Fig. 8: The ratio Eq. (44) for the u flavour, for X3 — ^3, as a function of ^3, at A^ = —0.15 
GeV^ (full hne), and at A^ = -0.25 GeV^ (dashed hne). 
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FIGURES 




FIG. 2. 
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FIG. 4. 
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